Improving quantum gate fidelities using optimized Euler angles 
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An explicit algorithm for calculating the optimized Euler angles for both qubit state transfer and 
gate engineering given two arbitary fixed Harriiltonians is presented. It is shown how the algorithm 
enables us to efficiently implement single qubit gates even if the control is severely restricted and 
the experimentally accessible Hamiltonians are far from orthogonal. It is further shown that using 
the optimized Euler angles can significantly improve the fidelity of quantum operations even for 
systems where the experimentally accessible Hamiltonians are nearly orthogonal. Unlike schemes 
such as composite pulses, the proposed scheme does not significantly increase the number of local 
operations or gate operation times. 
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I. INTRODUCTION 

Quantum computing [1] generally relies on the decom- 
position of arbitrary multi-qubit operations into products 
of elementary single and two-qubit gates, which must be 
implemented with very high fidelity. Although the avail- 
ability of an entangling two-qubit gate is crucial for uni- 
versal quantum computing [2], single qubit operations 
dominate virtually any decomposition of a multi-qubit 
quantum gate. For example, if we decompose a two- 
qubit gate into two-qubit gates that can be generated 
by a natural Ising interaction and local operations using 
the Cartan decomposition [3], at most three two-qubit 
terms are required in addition to 12 single qubit rota- 
tions. Therefore the fidelity of single qubit gates is crit- 
ical, as even small single qubit gate errors quickly accu- 
mulate, resulting in poor multi-qubit gate fidelities even 
if the entangling gate is perfect. 

One approach to improving gate fidelities and gate op- 
eration times is using optimal control. In general, opti- 
mal control fields can be derived by simultaneous opti- 
mization of many control parameters using numerical al- 
gorithms based on Poyntriagin's Maximum principle (see 
e.g. [4-9]). Optimal control may be the only viable option 
for implementing quantum gates for systems with highly 
complex Hamiltonians including off-resonant excitation 
and multi-body fixed coupling terms [10], but numerical 
optimization can be time-consuming and the resulting 
optimal control fields can be quite complicated and not 
necessarily easy to implement. By contrast, geometric 
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control, vaguely inspired by nuclear magnetic resonance 

(NMR) [11], requires only sequences of simple pulses to 
implement arbitrary single and multi qubit gates. Al- 
though compared optimally designed pulses the results 
may be suboptimal, this approach remains popular es- 
pecially in an experimental setting, due to its conceptual 
and experimental simplicity. However, there are limits to 
the applicability of standard techniques such as the Euler 
and Cartan decomposition, for instance, when we cannot 
implement local rotations about orthogonal axes, a situ- 
ation that arises in various settings, from global electron- 
spin achitectures [12] to charge-based semi-conductor 
quantum dot systems [13]. 

Geometric control generally relies heavily on Lie group 
decompositions such as the standard Euler decomposi- 
tion [14] of rotations in M^, which provides an explicit for- 
mula for decomposing any rotation in into a sequence 
of (at most) three rotations about two fixed, orthogo- 
nal axes, g and h. Due to the equivalence of S1[J(2) and 
SO(3) (SO(3) ~ §U(2)/{-l, 1}), it also provides an ex- 
plicit scheme to decompose any special unitary operator 
in S1U(2) into elementary complex rotations, combined 
with the generalized Cartan decomposition for multi- 
qubit gates, it provides a basis for generating arbitrary 
multi-qubit gates. The main drawback of the standard 
Euler angle decomposition is that requires orthogonal 
rotation axes, or respectively, Hamiltonians, while the 
Hamiltonians that are experimentally easily accessible 
are often at best approximately orthogonal, subject to 
certain simplifications such as negligible drift, rotating 
wave approximation, etc. Applying the standard Euler 
angle decomposition when the available basic Hamilto- 
nians are not orthogonal reduc;es the fidelity of most lo- 
cal gates, and hence virtually all multi-qubit quantum 
gates, regardless of the quality of the entangling gates, 
decoherence or other sources of noise that may reduce 
the fidelities of quantum operations. This is a significant 
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problem for applications such as quantum computation, 
where extremely high accuracy of the elementary gates 
is a prerequisite for scalability. 

One way to improve the accuracy of elementary gates 
is using composite pulse sequences [15 17] to compen- 
sate for certain systematic errors suc;li as rotation axis 
alignment and rotation angle errors. Such approaches 
have proved to be extremely valuable in ensemble-based 
quantum computing schemes such as liquid-state NMR, 
where multiple qubits are encoded into different nuclear 
spins of a larger molecule, and the system consists of 
an ensemble of a large number of identical molecules in 
solution. Due to magnetic field gradients, diffusion pro- 
cesses and inter-molecular interactions, the actual fields 
experienced by the individual molecules are subject to 
fluctuations, resulting in rotation angle errors, and to a 
lesser extent, rotation axis errors. Composite pulses re- 
duce these errors by replacing simple unitary operations 
(rotations) with sequences of rotations designed to "can- 
cel" certain errors. However, this systematic error can- 
cellation comes at the expense of increased overhead in 
the number of elementary operations, and hence time re- 
quired to implement a single quantum gate, especially if 
the systematic errors are so large as to require the use of 
concatenated composite pulses [18] . This can exacerbate 
other problems such as decoherence. Composite pulses 
can be applied to implement gates that are robust with 
respect to model uncertainty in non-ensemble-based sys- 
tems. However, unlike in ensemble-based schemes, where 
the systematic errors are a direct consequence of the fact 
that different molecules in the ensemble experience differ- 
ent forces, systematic errors due to model uncertainty in 
non-ensemble systems can be minimized by experimental 
system identification [19-22], and this has been shown to 
be advantageous in that it reduces to level of concatena- 
tion required for composite pulse sequences [18]. In this 
paper we show that if the actual Hamiltonians are known 
to sufficient accuracy then we can significantly improve 
gate fidelities with minimal overhead simply by optimiz- 
ing the Euler angles in the decomposition, potentially 
completely ehminating the need for expensive composite 
pulse sequences. 



II. QUANTUM GATE ENGINEERING USING 
LIE GROUP DECOMPOSITIONS 

Quantum computing generally relies on decomposing 
multi-qubit gates into products of elementary single and 
two-qubit gates, which can be applied simultaneously or 
sequentially to produce a desired unitary evolution. Fol- 
lowing the idea of using realistic physical Hamiltonians to 
generate quantum gates efficiently [23-25], one approach 
is to decompose a desired unitary operation U into ele- 
mentary unitary operations that can be easily generated 
by natural Hamiltonian flows of the system. For instance, 
given a system with a natural Ising coupling aPcrP be- 
tween adjacent qubits, and the ability to generate arbi- 



trary local unitary operations, any two-qubit gate can 
be factorized into a product of local operations, and the 
natural flows Z{t) = exp(— itai^Vi^') using the Cartan 
decomposition [26] 

U=Ui [Kl Z{ai)K,] [Kl Z{cx^) Ky\ Z{a^) U2, (1) 

where Z{a) corresponds to free evolution of the system 
under the Ising-coupling Hamiltonian for the time t = a, 
and Ui, U2, Kx and Ky are simultaneous local operations 
on both qubits. U\ and U2 depend on the particular gate 

to be implemented, while = U^\'k) ® u'f'\Ti), where 

U^\a) = exp(— i^ai'^^), fc = 1,2, and similarly for Ky. 
The Caxtan decomposition can be generalized to interac- 
tions involving more than two qubits [7], and an explicit 
algorithm to calculate the generalized Cartan decompo- 
sition was presented in [27]. Similar decompositions also 
exist for other natural non-local Hamiltonians but we 
still require very accurate single qubit gates. In principle 
such gates are easy to implement. Any W G §U(2) can 
be written as 

^ " sm(a)e ^ cos(q;) e '^y ^ ' 

where < a < f , < /3, 7 < 27r. We also have W = 
exp{—iH) with —iH e su(2), i.e., 

H = H{d) = dxCTx + dyay + d^a^ (3) 
with the usual Pauli matrices 

= (? J) ' 0') ' = (0 -1) ■ 

Let d = {d,,,,dy,d^) with U = ||d|| and n = fi'M. As 
= n^I, where I is the identity matrix, we have 

= cos(rit)I — i sin(rii)(n2,(Tj. -|- Uyay + nz<7z) 

_ f cos(r2t) — iuz sin(nf) — (% -|- iux) sin(r2f) \ 
~ \ (% ~ ^'"'x) sin(f2f) cos(fit) -|- iriz sin(f2f) J 

and comparison of the last equation with Eq. (2) shows 
that W = ex.p{—iTH) if we choose n, fl and T such that 

flT = arccos(cos a cos /3) (5a) 
n = — (sin a sin 7, sin a cos 7, cos a sin /3) (5b) 

with S = sin(rir). Thus, if we have full control over the 
single qubit Hamiltonians then we can implement any 
single qubit gate in a single step, and if there are no 
constraints on the magnitude fl of the Hamiltonian then 
the gate operation time T can be made arbitrarily small. 

Unfortunately, for most physical systems we cannot 
implement arbitrary Hamiltonians even locally. For ex- 
ample, the single qubit Hamiltonians for many potential 
qubit systems from ions to quantum dots are of the form 
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H = l{dzcr^_ + dx(Tx} OT H = \{dxax + dya^), restrict- 
ing us to rotations about axes in the xz or xy planes, 
respectively. If we have sufficient control over both 
and dz such as to be able to perform rotations about 
two orthogonal axes in the plane, then we can still im- 
plement arbitrary single qubit gates using the standard 
Euler decomposition, e.g., 

W{a, /?, 7) = [/, (/J + 7 - I) U,{a)Uz (/? - 7 + |) 

(6) 

where Ux{a) = exp(i^cra;) and Uz{a) = exp(i|^(T^) are 
elementary rotations about the x and z axis respectively. 

However, in practice there are often more constraints, 
limiting us to varying one or both parameters within a 
certain range. For example, for certain solid-state archi- 
tectures such as charge-based semi-conductor quantum 
dot systems [13], it is difficult or impossible to dynami- 
cally control the tunnel coupling d in the model Hamil- 
tonian H = Acuaz + da^- Thus d G [dmin-c?max] and if 
drnrn > then wc cannot implement rotations about the z 
axis, no matter how much control we have over the energy 
level splitting Aus, and practical constraints often make 
it impossible to find operating parameters (Al^-'i, di) and 
(Aw2,(^2) such that the corresponding Hamiltonians are 
exactly orthogonal. The same problem arises for other ar- 
chitectures where the amount of control is limited, such 
as global electron spin systems where many electron spins 
in quantum dots are simultaneously controlled by a fixed 
global field, and we can only control the detuning Aw 
of individual spins from the global field via local volt- 
age gates [12]. In other erases An..' may be fixed while 
we have limited control over the coupling strength or 
Rabi-frequency. 

In these examples (and other similar systems) we have 
a fixed drift Hamiltonian and constraints on a control- 
lable parameter. Without loss of generality, let us con- 
sider H{k) = I (da; + K<7z) with d > fixed and k e 
[0, Atmax]. Tr(a,a,) = and Tr{al) = Tr{al) = TV(I) = 2 
shows that the Hilbert-Schmidt inner product {H{ki) \ 
H{k2)) = T^[H{ki)^H{k2)] satisfies 

(P 
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{H{Kr) I H{k2)) < {H{0) I if(K„ax)) 

and \\H{k)\\ = y/{H{K) \ H{k)} = d^y{l + K^)/2. Thus, 
provided d ^ 0, the angle ( between the Hamiltonians 
H{0) and i?(Kmax) is determined by 



cos C = 



(g(0) I H,^. 



(8) 



Thus, C ~* f only for Kmax — * oo. For any finite value of 
Kmax the maximum angle between the accessible rotation 
axes will be less than |. If we use the standard Euler 
decomposition of a single qubit gate 



t/ = C/, (a) (7), 



(9) 



assuming U.j.(a) = cxp{—i'^X) and Uz{P) = exp(— i^Z), 
but the actual "z" -rotation is a rotation about Z = 



X cos ( + Zsin( with ( = (1 ± e) | then the gate actually 
implemented is 



U = U,{a)UI{l3)U,ij) 



(10) 



with C/|(/3) = exp(-zf Z). If there are no other errors 
the gate fidelity will be 

J-(/3,e) = ^|Tr(C/tc/)| = ^| Tr [£/,(/?) tc/,^ (/J) ]| ^^^^ 
= cos2(/3/2)[l -cos(e7r/2)] | cos(e7r/2)| 

and the gate error £(/?, e) = 1 — J-"(/3, e). Thus the max- 
imum single qubit gate error is 1 — | cos(e|)| = £{±Tr, e) 
|, the average single qubit gate 
cos(e^)], and for the maximum 
single qubit error to be below 10~^, the rotation angle 
error must be less than e = 0.9% or equivalently 



and, noting (cos^(a:;)) : 
error is favg(e) = hi^ 



cos 



sin^ 



2 

max 



> 0.9999 



(12) 



or Kmax > 70.7054. Hence, to keep the maximum gate 

error for a single qubit gate below the error threshold of 
10^^, for example, we would have to be able to make the 
energy splitting Aw (the controllable parameter) at least 
71 greater than the fixed coupling d, even if there were no 
other sources of error. If a CNOT-gate is implemented 
using the Cartan decomposition (1) with (Q;i,a2,a3) 



(7r/4,7r/4,0) and Ui = f/f^ (g) U^"' , C/; 



r(2) 



(1) 



u. 



(2) 



(1) 



(2) 



(1) 



Kj?\ where 







= ?7x(7r) 


(13a) 






= U^{-k)Uz{'k) 


(13b) 






= [/,(1.757r) 


(13c) 






= [/^(0.57r)C/^(1.57r)C/^(1.57r) 


(13d) 






= f74l.57r)f7^(0.57r) 


(13e) 



then assuming that our 2:-rotations Uz{/3) are really ro- 
tations about the tilted axis Z = Xsm{^e) + Zcos(|e), 
shows that the fidelity of the CNOT gate will be < 0.9999 
unless the rotation axis angle error is less than about 
0.6%, or Kniax ^ 100, cvcn if the entanglement-generating 
Ising-coupling terms are perfect and there are no other 
sources of error such as decoherence. In practice, other 
sources of error would mean that the error resulting from 
the rotation axis misalignment would have to be much 
smaller, and thus K,„ax much bigger, for the total errors 
to remain below the error threshold. Also note that for 
A^max = 1 the rotation axis angle error is 50%, and the 
maximum single qubit gate error is 1 — cos(7r/4), almost 
30%, and the error for a CNOT gate implemented using 
the Cartan decomposition above with unoptimized single 
qubit gates jumps to over 50%, assuming no errors in the 
Ising terms. 
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III. OPTIMIZED EULER DECOMPOSITION 

The previous section shows that accurate single qubit 

gates are crucial, and even small deviations of the ex- 
perimentally accessible single qubit Hamiltonians from 
orthogonality are problematic, not to mention situations 
where the experimentally accessible Hamiltonians are far 
from orthogonal. It is also known that any local unitary 
operation, in principle, can be generated exactly by per- 
forming a sequence of complex rotations about any two 
(fixed) Hamiltonians Hi and H2 that generate su(2), i.e., 
satisfy [Hi,H2] 7^ 0. Various Lie group decompositions 
have been considered for the related problem of imple- 
menting local qubit operations exactly in the presence of 
various types of fixed drift terms [28]. Ideally, however, 
we would like a simple explicit algorithm to calculate an 
optimal sequence of rotations given a fixed set of Hamil- 
tonians (rotation axes) and an arbitrary local gate. 

In the following we consider general decompositions of 
§0(3) instead of §1LJ(2) using the equivalence between 
SU(2) and §0(3) (modulo ±1). The advantage of con- 
sidering §0(3) is that it is easier to visualize rotations in 
R'^ than complex rotations in S1LJ(2). As a brief reminder 
we recall that any quantum state of a two-level system 
can be represented by a density operator 

P = p(r,^,</.)=2(^^2,.,3i^^ l-r^cos^j ^^^^ 

with 0<^<7r, 0<<;i<27r and < r < 1, and we can 
define a unique mapping between density operators p of 
constant purity Tr[p^] = i(l-t-r^) and points on a sphere 
of radius r in by 

(sin 6 cos (/)\ 
smOsincj) . (15) 
cos 6* ) 

The evolution of p under a constant Hamiltonian H then 
corresponds to a rotation of the Bloch vector s about the 
unit axis n = f2~^d with the constant angular velocity 
^ = ||d||, as defined in Sec. II. Given two Hamiltonians 
Hi and H2 we calculate the corresponding normalized 
Bloch vectors h and g and note that the angle between 
the rotation axes is given by 

C = arccos(h • g), (16) 

and C 7^ if and only if [iJi, i/a] ^ 0. 

Any target operator W{a,l3,"f) G §U(2) is equivalent 
(modulo ±1) to a rotation R{a,b,c) S SO(3) acting on 
the Bloch vector s with a ~ ^a, b = |(/3 + 7) and c = 
^(/J — 7), and we have explicitly 

(-61C1 + a2&2C2 &2Cl+a26lC2 -aiC2\ 
— 61C2 — a2&2Ci &2C2 — 026101 aici , 
0162 ai6i a2 J 

(17) 

where ai = sin(a), 02 = cos (a), and similarly for b and c. 



Algorithm 1 Calculate Generalized Euler Angles e for 
decomposition of arbitrary R G SO(3). 

Input: R G §0(3), unit vectors h,g G R^, h 7^ ±g. 

Output: Eulor angles e = (eo, . . . , 62^+2) 

C = arccos(h • g) 

(Oaf,<l>af) = Polar(i?h,h,g) 

p = \eaf/2Q 

if p > 

p = p-l 

end 

= eaf- Kp 

€1 = — arccos[cotCtan(0/2)] 

€2 = — arccos[(— cos^ C + cos ff) / ava? Q 

£3 = <t>af 

T, = (i?f,(7r)i?g(7r))^' 

Si = i?g(-ei)i?J-e2)Tpi?f,(-e3) 

(fbi, 061) =PoleLr{SiRRy{n/2)h, h, g) 

€0 = 4>bl 

if p = 3 

e = [eo, ei, £2 -|- £3] 

else 

e = [eo, ei, e2 -h 71, 71, . . . , 71, es] 

2p-l 

end 

e = e-|-27iL(27i-e)/(27i)J 



Algorithm 2 Polar: Compute polar coordinates of a 
with respect to orthogonal frame induced by h and g 

Input: unit vectors a, h, g G K"^ 
Output: polar coordinates {6, 4>) of a 
z = h 

y = h X g 

y = y/l|y|| 

X = y X h 

9 = arccos(a ■ z) 

(j) = arctan(a ■ y, a • x) 



Proposition. Any rotation R e SO(3) about an arbi- 
trary axis in can be decomposed in a series of rota- 
tions about two (non-identical) fixed rotation axes h and 
g in as follows: 

o_ / ^h(e3 + e2)i?g(ei)i?hM, 

\ Rf^{es)XP-'R^{n)Ri^{i2)R^{ei)Ri^{eo), p > ' 

(18) 

where the parameters p, eg, ei, £2 and e^, are given explic- 
itly in Algorithm 1, and X — R^{'k)R-^{it), £2 = £2 + tt. 

We would like emphasize here that the important part 
from an application point of view is not the existence 
of a general decomposition of the form (18), which was 
predicted by [29] and shown in [30] , but the simple Algo- 
rithm 1 to compute the Euler angles e = (eo, . . . , £2^4-2) 
in the decomposition based on analytical formulas. In 
most cases the Euler angles in decompositions about non- 
orthogonal rotation axes can only be determined numer- 
ically using optimization techniques [31], but in this case 
we are in the fortunate position that we can derive rela- 
tively simple analytic formulas for all Euler angles. 
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Before wc discuss the derivation of this result and the 
algorithm, we should briefly justify the use of the expres- 
sion "optimized" Euler angles in the title. The factoriza- 
tion (18) shows that in general 2(p — 1) + 5 = 2p + 3 rota- 
tions are necessary, and for ^/ < tt the maximum number 
of steps is 2p + 3 < [|^] + 1, which is equal to the order 
of generation of §0(3), which is A; = \^] + 1 according 
to Lowenthal's criterion [29]. This means that the de- 
composition is optimal in the sense that we achieve unit 
fidelity, and that we cannot generate arbitrary rotations 
using rotations about the fixed axes h and g in fewer 
steps in general [34]. Optimality in terms of the num- 
ber of rotation steps is often related to time-optimality 
as more rotation steps generally will take longer to com- 
plete, although time-optimality is not guaranteed. For 
instance, if rotations about one axis can be implemented 
much faster than rotations about the other then a se- 
quence that requires more steps but fewer slow rotations 
by larger angles may be faster to implement. Also, it 
should be noted that while the decomposition (18) gen- 
erally provides the best way to implement a quantum 
gate in S1[J(2) exactly using a minimal number of rota- 
tions about two fixed axes, we may be able to implement 
a particular gate substantially faster if we can dynami- 
cally vary the rotation axes continuously. For example, 
consider a system with Hamiltonian H = ^{(Tx + kGz)- 
If we can temporally vary n{t) to take any value in the 
range [0, Kmax] rather than two fixed values, e.g., and 
'^max , then we may be able to implement a particular gate 
faster by numerically optimizing K,{t) instead of using the 
optimized Euler angle decomposition. The attractivity of 
the generalized Euler decomposition lies in its simplicity. 
It is a simple "bang-bang" control scheme that can be 
used to achieve unit fidelity when we have limited control 
and cannot (or do not wish to) implement complicated 
temporal control field profiles K{t). 

The derivation of the algorithm is based on steering of 
a state represented by a vector sq e of length r to 

another state sj G K"* the same distance from the origin 
using only rotations about the two fixed rotation axes 
given by the unit vectors h and g. An algorithm and 
detailed explanation how to move from one point on the 
sphere to another though a sequence of rotations about 
two fixed axes h and g is presented in Appendix A. How- 
ever, this algorithm on its own is not sufficient for gate 
engineering as a single point and its image on the unit 
sphere in are not suflticient to uniquely determine a 
rotation in i? G SO(3). Rather, we need at least two 
(non-antipodal) points and their images to fix R. This 
may seem surprising as the parametrization (2) for a uni- 
tary operator W e S1[J(2) shows that the image W\'tp) of a 
single Hilbert space vector ji/') is sufiicient to fix all three 
parameters. The corresponding real rotation (17), how- 
ever, cannot be fully determined by the image of a single 
Bloch vector. Rather the mapping of a single point in 
M.^ defines a one-parameter family of elements of SO(3), 
and a second (non-antipodal) point and its image are re- 



quired to fix all three parameters a, 6, c of R. The reason 
for this apparent discrepancy is that when we transform 
from complex Hilbert space vectors \^) to density matri- 
ces p = the information about the global phase of 
the initial and final states, which helps fix W uniquely, 
is lost. 

Any pair of non-antipodal initial points (ao,lJo) and 
their images (a/,lj/), are sufiicient, but it is convenient 
to choose the initial points to be the orthogonal set of 
initial states 

ao = h=(0,0), (19a) 
bo = i?y(f)ao= (f,0), (19b) 

where the pairs {9, (j)) are the polar coordinates 

6 = 9{a) = arccos(a • z) (20a) 
^ = ^{a) = arctan(a • y, a • x) (20b) 

of the unit vectors ao and bg with respect to the rectan- 
gular coordinate system (x, y,z) defined by 



y = 



h X g 



|h X g|| 

and the four-quadrand arctangent is 



y X h. 



(21) 



arctan(i/, x) = < 



' arctan \y/x\ x > 0, y > 

TT — arctan \y/x\ x < 0, y > 

n + arctan \y/x\ x < 0, y < 

, 27r — arctan \y/x\ x > 0, y < 0. 



(22) 



Then we compute the polar coordinates of the corre- 
sponding images under the target rotation R 



a/ 



i?ao = 



{9af,(f>af), 



(23a) 
(23b) 



and use the state transfer algorithm (Algorithm 3) to 
calculate Si, the series of rotations that steer aj to ao, 

5i = i?g(-ei)i?j^(-e2)rPi?h(-'^a/) 

where T = R^{TT)Rg{TT) and p is the largest integer 
strictly less than 6'a//2C, i.e., p = lOaf/^C] - 1. The 
generalized Euler angles ci and £2 can be obtained by 
inserting Of = and 6*0 = Oaf — 2pC 7^ into subroutine 
PR2 (see Appendix B) 



€2 = — arccos 



ei = — arccos 



2pC 



cot C tan 



- COS^ C + COsjOaf - 2pC) 

sin^C 



(24a) 
(24b) 



Since Sihf = (^,(f>hi) and bo are unit vectors with 
the same 0-angle {0 = |), the same series of rotations 
preserves the distance between the points, hence 

bo = i?h(-0bi)^ib/. 
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The rotations leave do unchanged and we have 

ao = Rf^{-(l)bi)Sia.f. 
Thus R = /S'f^i?j^(^bi) and we have the decomposition 

R = i?f,(e3)Xf i?j^(e2)i?t(ei)i?h(eo) (25) 

where X = T'^ = i?g(7r)i?^(7r) and eo = (j)bi, es = 4>af 
and ei, as in (24). If p > 1 we can combine the two 
subsequent h rotations, while for p = 0, = I, thus the 
optimal decomposition is 

o _ / ^h(e3 + e2)i?t(ei)i?h(eo), P = 

\ R^{e,)XP-'R^{7T)R^{h)R^{ei)Ri,{eo), p > 

with €-2 = €2 + TT, which completes the proof. 



gate increases from w 10""* for k = 100 to over 51% for 
K = 1. Also note that the penalty for non-orthogonal 
Hamiltonians in terms of the number of rotation steps 
required is actually rather small unless k is very small. 
Indeed for k > 1 all of the elementary gates in the table 
can be implemented in at most four steps, and for k > 1, 
this is indeed the maximum number of steps required for 
any single qubit gate. To see this recall that Lowenthal's 
criterion guarantees that the maximum number of steps 
in the decomposition of any single qubit gates is 



K = 









+ 1 = 


c 





arccos[(l 



-1/21 



+ 1, (30) 



which yields K = 3 ioi k = oo and K — A for 1 < k < oo, 
iV- = 5 for K = 1, /ST = 6 for ^[cos(7r/5)]-2 - 1 < k < 1, 
and so forth. 



IV. APPLICATIONS 

To apply the results in the previous section to imple- 
ment a quantum gate 

U = exjp[i^{nxax + Uydy + n^cFz)] (26) 

given the Hamiltonians 

Hi = ^a^, H2 = ^{a^ + na^), (27) 

we identify the normalized Hamiltonians Hi = <Jx and 
H2 = {cTx + k<Tz)/vT+~k2 with the unit vectors h = 
(1,0,0)-^ and g = (1, 0, k)-^/\/1 + k^, respectively and 
use Algorithm 1 to decompose the corresponding SO(3)- 
representation of the target operator U 

A = exp[$(n^i?^ + TiyRy + n^R^)] (28) 

where Rx, Ry and R^ are the rotation generators 

/O 0\ /O -2\ / 2 0\ 

Rx=\^ 2 ,i?y= 00 ,i?^= -2 00 
\0-2 0/ \2 00/ \000/ 

Table I shows the optimized Euler angle decomposition 
results for the gates S = exp (i^tTx), T = exp (if fx), 

f/Had = exp (i^(a, + a,)) (29) 

as well as for the single qubit gates (13) required to im- 
plement a CNOT gate via the Cartan decomposition (1). 
Gates that require only ctx rotations have been omitted 
as they are trivial to implement with the given Hamil- 
tonians. Using the optimized Euler angles rounded to 
four significant digits, the gate errors for all single qubit 
gates in the table, as well as the CNOT gate, are be- 
low 3 X 10 ~^ for values of k shown, while the gate errors 
using the standard Euler angles increase to almost 30% 
for Ky and k = 1. The error for the resulting CNOT 



V. CONCLUSIONS 

The Euler decomposition of unitary operators in S1[J(2) 
is widely used to implement single qubit gates by de- 
composing them into products of rotations about two 
orthogonal axes determined by fixed Hamiltonians. The 
approach can be problematic however as experimentally 
accessible Hamiltonians in many cases may not be or- 
thogonal. Depending on the situation, in some cases the 
Hamiltonians can be made almost orthogonal, while in 
others the constraints may be far more severe. In either 
case, however, lack of orthogonality of the underlying 
Hamiltonians leads to errors in the gates implemented, 
and even small errors can propagate. A rotation axis an- 
gle error of even 1% results in single qubit gate errors 
above the error threshold of 10~^ even if there axe no 
other sources of errors, and the single qubit errors com- 
pound and lead to even larger errors for two-qubit gates. 
Such systematic errors can easily be corrected, however, 
by adapting the Euler decomposition to the actual Hamil- 
tonians available. 

We have presented a explicit algorithm to calculate the 
optimized Euler angles for any single qubit gate and two 
arbitrary fixed Hamiltonians, and shown that we can sub- 
stantially improve single and two-qubit gate fidelities by 
using optimized rather than standard Euler angles. The 
idea is attractive because the computational overhead to 
calculate the optimized Euler angles is minimal and the 
implementation is no more demanding than standard ge- 
ometric control, i.e., no additional resources are required. 
There is a small price to pay in terms of an increase in 
the number of rotation steps required to implement a par- 
ticular gate, but unless the maximum angle between the 
experimentally accessible Hamiltonians is very small, this 
increase is very slight, e.g., from at most three steps for 
orthogonal Hamiltonians to four for Hamiltonians with 
angle ( greater than 45° and five if C = 45°. For a model 
Hamiltonian H{k) = ^{(Tx + i^'^z) with a fixed coupling 
parameter d, this condition is satisfied if the energy level 
splitting can be made at least as large as the tunnelling 
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K 


f„ (%) 




<H2) 






T 


oo 








1.7500 









100 


0.0007 


0.0013 


1.7500 


0.0013 






50 


0.0029 


0.0026 


1.7499 


0.0026 






10 


0.0727 


0.0132 


1.7487 


0.0132 






5 


0.2844 


0.0264 


1.7448 


0.0264 






1 


4.2893 


0.1359 


1.6359 


0.1359 




S 


oo 








1.5000 









100 


0.0025 


0.0032 


1.5000 


0.0032 






50 


0.0100 


0.0064 


1.4999 


0.0064 






10 


0.2481 


0.0319 


1.4968 


0.0319 






5 


0.9710 


0.0641 


1.4873 


0.0641 






1 


14.6446 


0.5000 


1 


0.5000 




t^Had 


oo 





1.5000 


1.5000 


1.5000 






100 


0.0025 


1.5032 


1.5000 


1.5032 






50 


0.0100 


1.5064 


1.4999 


1.5064 






10 


0.2481 


1.5319 


1.4968 


1.5319 






5 


0.9709 


1.5641 


1.4873 


1.5641 






1 


14.6442 





1 









oo 





0.5000 


1.5000 


1.5000 






100 


0.0025 


0.5032 


1.5000 


1.5032 






50 


0.0100 


0.5064 


1.4999 


1.5064 






10 


0.2481 


0.5319 


1.4968 


1.5319 






5 


0.9709 


0.5641 


1.4873 


1.5641 






1 


14.6443 


1 


1 









oo 








1.5000 


0.5000 






100 


0.0025 


0.0032 


1.5000 


0.5032 






50 


0.0100 


0.0064 


1.4999 


0.5064 






10 


0.2481 


0.0319 


1.4968 


0.5319 






5 


0.9709 


0.0641 


1.4873 


0.5641 






1 


14.6445 


0.5000 


1 


1 






oo 











1 


1 




100 


0.0050 


0.5000 


1.9936 


0.5000 


1 




50 


0.0200 


0.5001 


1.9873 


0.5001 


1 




10 


0.4963 


0.5032 


1.9362 


0.5032 


1 




5 


1.9419 


0.5127 


1.8718 


0.5127 


1 




1 


29.2893 


1 


1 


1 


1 



TABLE I: Optimized Euler angles (in units of tt) for various 
single qubit gates and different values of k. e(-ffi) indicates 
a rotation by e about the normalized axis Hi. So is the gate 
error that results if the standard Euler angles for k = oo 
are used. The gate errors using the optimized Euler angles 
truncated to four decimal digits are < 3 x 10^® for all gates, 
and can be made arbitrarily small by increasing the number 
of significant digits of the Euler angles. 



energy d, or k = 1, whereas the standard Euler decom- 
position would require k ^ oo, or energy level splittings 
that are orders of magnitude greater than the tunnelling 
energy d to achieve near-orthogonal Hamiltonians. 

Overall, the overhead in terms of complexity of the 
pulse sequences is small compared to alternative ways to 



correct for rotation axis errors, such as composite pulse 
sequences, and this overhead seems acceptable, consid- 
ering that relaxing the need to be able to perform ro- 
tations about orthogonal axes may allow for substantial 
simplifications of the underlying architectures. Another 
source of overhead of the technique is the need for ini- 
tial characterization of the Hamiltonians. It must also be 
stressed that optimized Euler angles are designed to min- 
imize errors for a single system. They cannot compen- 
sate for ensemble errors, i.e., errors arising from the fact 
that individual systems in a large ensemble may experi- 
ence different fields and thus different effective rotations. 
However, the approach is an effective way to improve 
gate fidelities for non-ensemble systems with non-ideal 
Hamiltonians. 

Further work is necessary to extend the results to 
higher-dimensional systems. Another issue is that differ- 
ent gates require different amounts of time to implement. 
This is not a problem for a single system but would be 
for a large register if one wants to implement gates on 
different qubits simultaneously. Here we have only used 
two fixed Hamiltonians with a fixed angle between them. 
In many cases, however, we may be able to vary the con- 
trollable parameter continuously up to some maximum 
value. An interesting question in this regard is whether 
we can exploit the (limited) variation in the tilt angle to 
design simple geometric controls that allow us to imple- 
ment arbitrary gates in a fixed amount of time. 
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APPENDIX A: STATE TRANSFER ALGORITHM 

The objective of state transfer is to steer the system 
from a known initial state Sq — {r, Oo,(f)o) on a sphere of 
radius r to a target state Sf = {r' ,9f,(j)f). With uni- 
tary control only states on the same sphere as the initial 
state are accessible by performing a sequence of rotations 
about the axes h and g, respectively. We shall assume 
r = 1, noting that the sequence of rotations that steers 
the normalized initial state ag to the normalized target 
state SLf steers sq to s/ if they lie on the same sphere of 
radius r, and to a state that is as close to the target 
state as we can get with unitary control if r ^ r'. 

A general strategy to get from Sq to Si with a min- 
imum number of rotations about the axes h = z and 
g = (sin ^,0, cos C)-^, following [30] and earher work [29, 
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FIG. 1: Sphere with arbitrary rotation axes, respective coor- 
dinate system and an arbitrary vector a (angles 6 and (j))- 



32, 33], is to rotate the initial state by a suitable angle 
about either axis to map it to a point on the great circle 
in the x, z plane, followed by a sequence of 7r-rotations, 
alternating about the h and g axis, until the angle 0' of 
the current state differs by less than 2C, from the Of oi 
the target state, i.e., we are within direct reach of the 
target state, followed by a final rotation by a suitable 
angle about the same axis we started with. By Lowen- 
thal's criterion any state can be reached from any other 
state in at most fc + 1 steps, k being the smallest inte- 
ger > ^ — 1 [29], and for C = f at most two steps are 
required. 

Based on this idea, we can derive an explicit algorithm 
for calculating the generalized Euler angles of an optimal 
decomposition given the angle a between the rotation 
axes h and g, and the relative coordinates (0O:0o) ^^id 
{9f,(t)f) of the initial and final state. Assume 6*0 > Of 
and C < 7r/2. If Oq — Oj < 2( then we can get from {Oq, 0) 
to {Of,0) in two steps, either by rotating (6*0,0) around 
g by an angle 0, followed by a rotation about h by an 
angle (j) (Subroutine PRl), or by a rotation around h by 
an angle (p, followed by a rotation around g by an angle 
(Subroutine PR2) . li Og — Of > 2( then we move from the 
initial point to a point with Oq — Of < via a sequence 
of TT-rotations around axes g and h as described before. 
If C < tt/2 but Oq < Of then we exchange the initial 
and final points, apply the algorithm and finally reverse 
the sequence of rotations. If C > we set (6*1, = 
(tt - 6*0, (6*2, (1)2) = (tt - Of, cj)f) and C = tt - ( and 
apply the algorithm. 

The algorithm returns a list of pairs (e,r), where e e 
[0,7r] is a generalized Euler angle and r = h or r = g 
indicates the rotation axis, which defines the necessary 
sequence of the rotations. E.g. if 6*0 > 0/, C < f f^^d 



routine PRl was used, then 

s{Of,d^f) = Ri^{c^f)Ri^{mmT''Ri,{-M<OQ,cpQ). 

(Al) 

where the angles <f) and are given by 

sin Oq — cot C(cos r — cos Oq) , , ^ s 

cos (j) = — (A2a) 



cosO 



sm u f 

- cos CC0S(C — Oq) + COS Of 



sinCsin(C — Oq) 
Similarly, if PR2 was used 



(A2b) 



s{Of,^f) = i?f^(</./)i?g(0)i?j^((/))TPi?^(-~0o)s(0o>o) 

(A3) 

where the angles (j) and are determined by 

sinflf — cotCfcos^o — cos^f) , 

cos = — — (A4a) 

sm 6^0 

cos = -^°^CcoB(C-^/)+cosgo 

sinCsin(C - 6*/) ^ ' 

The procedures PRl and PR2 are described in Ap- 
pendix B, and in both cases we have 

TP = Ri,{n)R^{TT) . . . Ri^{7r)R^{n) . (A5) 
2p rotations 

The minimal number of steps can be calculated explic- 
itly. If 6*0 = 6*/ then only a single rotation about h is 
required, otherwise the minimal number of steps to get 
from (Oq, 0o) to (Of, (pf) starting with a rotation about h 
is iV = 2p -t- 2 + Q{q), where Q{q) = 1 for g > and 
otherwise, 

p - int[(0o - 0f)/2C] (A6a) 
q~cosOf~ [tan C(sin — sin6'/ cos 0/) -I- cos ^] (A6b) 

and ~ Oq ~ 2pC,. Here, int(x) indicates the integer part 
of X. 

The minimal number of steps to get from (Oo,(f)o) to 
{Of,(j)f) starting with a rotation about g is iV' = 2 -|- 
6(61' - Of) [2p+l + e{q)], where 

p = int[(0' - 6l/)/2C] (A7a) 
q — cosOf~ [tan C(sin — sin6'y cos 0/) -t- cos ^] (A7b) 

with 0' — arccos(cos C cos Oq + sin ( sin Oq cos 4>o) — C and 
0^0'- 2K- 



APPENDIX B: SUBROUTINES PRl AND PR2 

Both procedures take only the 0-angles of the initial 
and final points, Oq and Of, respectively, and the angle 
a between the axes h and v as input, assuming that the 
points have already been shifted to the x — z-plane. 
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Algorithm 3 Calculate sequence of rotations about 
arbitrary, fixed axes h and g required to move from one 
point on the unit sphere to another. 

Input: {80, <j>o), {6f,(j)f) - polar coordinates of initial and 
final point with respect to relative coordinate system (21). 
Output: List of pairs (e,u), u £ {h, g} defining rotation 
steps necessary to get from initial state to final state, 
if 6^0 ~ Of return {(0/ — (/)o,h)} 
else 

p= L(^o-e/)/2CJ 

80^60- 2pC 

if 6*0 = ^/ return 

{(<^/, h), (tt, g), (tt, h), . . . (tt, g), (tt, h), (-(jio, h)} 

^ V ' 

2p pairs 

else if 00 < C or (2C - 6*0 > and 2^ - > 9f) 
return 

2p pairs 
PR2[0o,e/,C],(-0o,h)} 

else return 

{(<^/. h), (tt, g), (tt, h), . . . , (tt, g), (tt, h), 
2p pairs 



end 



end 



end 



Subroutine PRl 



Fig. 2 shows that ||MB|| = ||MQ'|| and ||NA| 
IINQ'II and 



IMB-MQ'I ^ INA-NQ'I 

cos = 





FIG. 2; Rotations on the sphere and projection to the x — z 
plane for subroutine PRl 



IMBI 



INAI 



(Bl) 



Noting that A — (sin 6*0, 0, cos 0o) and Q' = {q,p,cos9f) 
for suitable values of p and g, and taking Q = 
(g, 0, cos 0/) to be the projection of Q' onto the xz plane, 
shows 



cot ( = — 



sm t/o 



QA i 

Q A • z cos ^0 — cos I 



(B2) 



cos(0o — C) and |NA| = sin(0o ~ Oi and 

NA = (sin 00 — r sin C, 0, cos Oq — r cos C) 
NQ' = {q — r sin (, p, cos Oq ~ r cos C) 
NA • NQ' = (sin 6*0 - n-x){q - n^) + (cos 6*0 - n,f, 



and thus 5 = sin 6*0 — cot C (cos df — cos6o). Noting further 
that B = (sin , 0, cos 6*^) and M = (0, 0, cos 6'y), shows 
that MQ' = {q,p,0), MB = (sin6'/, 0, 0), and therefore 



which after some simplification gives 



COS0 = 



qsinOf sin 6*9 — cot C( cos 61/ — cos 6*0 ) 



(B3) 



cos e = -^oKcos(C^go)+cosg, 
sin ^ sin (( — Oq) 



Furthermore, we have N = r(sin C, 0, cos C) with r = 
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FIG. 3; Rotations on the sphere and projection to the x — z 
plane for subroutine PR2 



Fig. 3 shows that ||MA 
INQ'll and 



Subroutine PR2 

IIMQ'II and ||NB| 



cos ( 



|MA||MQ'| 
|MA|2 



cos f = 



|NB||NQ'| 
INBP 



(B5) 



Noting that B = (sin 6*^ , 0, cos 6*^) and Q' = ((7,p, cos 6'o) 
for suitable p and q as before, shows that 



cot C = — 



QB i 
QB z 



sin 6f — q 
cos 9f — cos 9o 



(B6) 



i.e., q = sin0/ — cotC(cos0o — cos 6'/). 

Taking Q — (g, O,cos0o) to be the projection of Q' 
onto the xz plane, and noting that A — (sin 610, 0, cos 6*0) 
and M = (0,0, cos 6*0) shows that MQ' = ((?,P,0), 
MA (sine*/, 0,0), and thus 

gsin^Q sin^f — cotC(cos0o — cos6'f) , „, 

cos = „ = i ^. , (B7) 

sin 9o sm 6^0 

Furthermore, we have N = r(sin C, 0, cos C) with r = 
cos(6l/ - C) and |NB| = sin(6'/ - C), and 

NB = (sin 6*/ — r sin (, 0, cos df ~ r cos C,) 
NQ' = [q — r sin p, cos 9q ~ r cos C) 
NA -NQ' = (sin 6*0 -n,^){q-n^) + {cos9a - n^^f , 

which after some simplification gives 



cos 9 — — cot C cot (C — 9f) + cos 00 cscCcsc(C — 9f) 
— cosCcos((^ — 9f) + cos6'o 



sin ('sin(C — 9f) 



(B8) 
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